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We suggest a mechanism that may resolve a conflict between the precession of a neutron star and 
the widely accepted idea that protons in the bulk of the neutron star form a type-II superconductor. 
We will show that if there is a persistent, non-dissipating current running along the magnetic flux 
tubes the force between magnetic flux tubes may be attractive, resulting in a type-I, rather than 
a type-II, superconductor. If this is the case, the conflict between the observed precession and 
the canonical estimation of the Landau-Ginzburg parameter k > l/v2 (which suggests type-II 
behaviour) will automatically be resolved. We calculate the interaction between two vortices, each 
carrying a current j, and demonstrate that when j > where q is the charge of the Cooper pair 
and A is the Meissner penetration depth, a superconductor is always type-I, even when the cannonical 
Landau-Ginzburg parameter k indicates type-II behaviour. If this condition is met, the magnetic 
field is completely expelled from the superconducting regions of the neutron star. This leads to the 
formation of the so called intermediate state, where alternating domains of superconducting matter 
and normal matter coexist. We further argue that even when the induced current is small j < 
the vortex Abrikosov lattice will nevertheless be destroyed due to the helical instability studied 
previously in many condensed matter systems. This would also resolve the apparent contradiction 
with the precession of the neutron stars. We also discuss some instances where anomalous induced 
currents may play a crucial role, such as in neutron star kicks, pulsar glitches, the toroidal magnetic 
field and the magnetic helicity. 

PACS numbers: 97.60.Jd, 26.60.+C, 74.25. Qt, 97.60.Gb 



I. INTRODUCTION 

This paper is motivated by calculations [l[ that show 
a conflict between the observed magnitude (~ 3°) and 
frequency (~ 1 per year) of neutron star's precession [2| 
and the widely accepted idea that protons in the bulk 
of the neutron star form a type-II superconductor. We 
begin with a review of this contradiction. 



A. Precession and Superconductivity 

In the generally accepted picture, the interior of a neu- 
tron star contains neutrons and a small number of pro- 
tons and electrons. A compilation of neutron and proton 
scattering data [3| implies that the extremely cold (10 s 
K), dense (10 15 g/cm 3 ) nature of the neutron star should 
cause the neutrons to form 3 P 2 Cooper pairs and become 
a superfiuid and the protons to form 1 Sq Cooper pairs 
and become a superconductor. 

A rotating superfiuid cannot form a solid body to carry 
circulation but instead forms vortices of quantized circu- 
lation that run in the direction of the angular velocity 
It is also generally accepted that the protons form 
a type-II superconductor, which means that it supports 
a stable lattice of magnetic flux tubes in the presence of 
a magnetic field ||. The rotation of the neutron star 
(1 — 10 3 Hz) and the presence of the enormous magnetic 
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field (~ 10 12 G) are sufficient for both superfiuid and su- 
perconducting vortices to form inside the neutron star. 
It should be noted that when a superconductor is placed 
in a rotating container it co-rotates with the container 
at the expense of a small current known as the London 
current [6j. 

Vortices formed in a superfiuid will always repel and 
form a lattice. Vortices formed in a superconductor can 
either attract, in which case the magnetic field is expelled 
from a superconductor, or they can repel each other and 
form a triangular lattice. These two behaviours label the 
superconductor. If the vortices attract, the superconduc- 
tor is called type-I and if they repel it is called type-II. 
The distinction between the two types of superconduc- 
tivity will be very important in this paper. 

The standard way to determine the type of super- 
conductivity is to calculate the value of the Landau- 
Ginzburg parameter k = j, where A is the Meissner 
penetration depth and £ is the coherence length of the 
superconductor. If k > l/\/2 then the superconductor is 
type-II, otherwise it is type-I. This parameter is weakly 
dependent on density. Typically for nuclear matter in 
a neutron star, A ~ 80 fm and £ ~ 30 fm [![, which 
means that K ~ 2.6. At higher densities it is possible to 
have a Landau-Ginzburg parameter that indicates type-I 
behaviour [7[, but only in the densest part of the core. 
There still exist large regions of the neutron star where 
type-II behaviour is predicted. 

In a neutron star both condensates are subjected to 
angular momentum and magnetic flux and form lattices, 
but the proton vortices (~ 10 24 per m 2 ) are much more 
numerous than the neutron vortices (~ 10 10 per m 2 ) and 
are tangled around them. It is the formation of these 
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lattices that cause the contradiction with the precession 
of the neutron star. 

The existence of precession means that the superfluid 
neutron vortices no longer form along the rotational axis 
of the star, but along the axis that is the sum of the pre- 
cession and angular momentum vectors. When the star 
precesses the vortices now move with respect to the rota- 
tion of the star and, in turn, with respect to the proton 
vortices which are entangling them. If the precession is 
large enough one of two things must happen; either the 
neutron vortices move with the proton vortices or they 
pass through each other. 

Requiring the neutron and proton vortices to move to- 
gether places severe restrictions on the precession. This 
is the case of "perfect pinning" that was discussed in 
using macroscopic dynamics. We will follow the argu- 
ments in [l[ as we are specifically interested in proton 
vortices as a mechanism for pinning. Because the core of 
the star is superconducting, the proton vortices, which 
carry magnetic flux, are resistant to being moved Q and 
thus the neutron vortices are restricted to move slowly. 
This means that the neutron vortices are pinned to the 
rotation of the protons and thus are pinned to the rota- 
tion of the crust. If this pinning is present the neutron 
star can only precess at very high frequencies. If the star 
is to precess more slowly at large amplitudes then it is 
necessary for the neutron vortices to pass through the 
proton vortices. 

The case of " imperfect pinning" was first discussed in 
fTol ] using the concept of vortex drag. Microscopically, 
this drag is created by large numbers of neutron vor- 
tices passing through proton vortices [l[. This creates a 
number of excitations and is a highly dissipative process. 
Both methods find that the precession is highly damped 
and that there are no modes of large, persistant preces- 
sion. Based on these estimates it is concluded that given 
the observed precession neutron vortices and proton flux 
tubes cannot coexist in the star Either the star's 
magnetic field does not penetrate any part of the core 
that is a type-II superconductor or that at least one of 
the hadronic fluids is not superfluid. Based on pairing 
calculations that predict neutron and proton superfluids 
coexist in the outer core the latter is very unlikely, so we 
will look to the former. 

If the core is a type-I superconductor, the magnetic 
flux could exist in macroscopic regions of normal mat- 
ter that surround superconducting regions known as an 
intermediate state. The magnetic flux would not form 
proton vortices and there would be nothing to impede 
the movement of the neutron vortices, thus allowing the 
star to precess with a long period. 

For proton vortices to not form a lattice, a mecha- 
nism must be present to make the interaction between 
them attractive (type-I behaviour). If it does so even 
when the cannonical Landau- Ginzburg parameter sug- 
gests that the superconductor is type-II then the incon- 
sistency described above will be resolved. 



B. A New Mechanism for Type-I 
Superconductivity 

In this work we suggest a mechanism that leaves the 
Landau-Ginzburg parameter unchanged, but causes the 
system to behave quite differently than in the standard 
picture. To be precise, we will show that even when 
k > the system prefers the intermediate state and 
that the apparent contradiction between n ~ 2.6 and the 
observed precession is avoided. 

This is achieved if the system supports a persistent, 
non-dissipating current running along the core of a vor- 
tex. Such topological currents appear in many systems as 
the consequence of a quantum anomaly. It is well known 
that anomalies, and the topological currents they induce, 
have important and non-trivial implications; the electro- 
magnetic decay of neutral pions 7r° — ► 27 is a textbook 
example. 

As we will discuss later, analogous topological currents 
have even been observed in some condensed matter sys- 
tems. However, with a few exceptions, the analysis of 
quantum anomalies has not received attention in the lit- 
erature devoted to dense matter systems in general, and 
neutron stars in particular. In section II we discuss in 
detail why and when such topological currents may arise 
in high density systems. 

If such currents are induced, they drastically change 
the behaviour of the system. Normally the interaction 
between superconducting vortices has two terms: an at- 
tractive force which comes from the order field and a 
repulsive force which comes from the gauge field. In the 
presence of an induced current, a third, attractive force 
will appear and, if the current has sufficient magnitude, 
the system will behave as a type-I superconductor. This 
situation will be discussed in sections II and IV. 



C. Relation to Previous Work 

The same problem has been discussed previously in 
[Tl| , [l2j . In particular, has shown that the existence 
of a type-I superconductor in a neutron star resolves the 
conflict addressed in Reference [III starts by assum- 
ing that the equilibrium structure of a type-I supercon- 
ductor contains alternating superconducting and normal 
domains. By use of a hydrodynamic restriction based 
on the moment of inertia of the crust and the moment 
of inertia of the superfluid it is shown that the alternat- 
ing domain structure seen in type-I superconductors will 
always allow for undamped precession. 

This result can be understood in a different way. The 
arguments of reference [l| relied on the proton vortices 
tangling around the neutron vortices. In a domain struc- 
ture there is more room for the neutron vortices to move 
unhindered by the proton vortices allowing for large am- 
plitude, high frequency precession. 

However, the calculation in reference [ll| starts by as- 
suming that type-I superconductivity already exists in 
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neutron stars. It does not investigate how a type-I super- 
conductor could arise and it is not obvious how to recon- 
cile this with standard arguments that suggest k ~ 2.6, 
which ambiguously implies that the superconductor is 
type-II in a finite volume of the neutron star. 

In reference [l2| a mechanism has been suggested that 
potentially resolves the conflict between the standard es- 
timation of k and type-I behaviour. The mechanism is 
based on the idea that proton and neutron Cooper pairs 
have almost identical interactions, akin to their underly- 
ing isotopical symmetry, even though the Fermi momenta 
and densities for protons /z p and neutrons ^ n are vastly 
different. A small difference in interactions was modeled 
by a small effective asymmetry parameter £ <C 1 • In this 
case it has been shown that the type of superconductivity 
is not governed by the canonical Landau- Ginzburg pa- 
rameter k, but instead, by an effective Landau-Ginzburg 
parameter k = £k. This effective parameter can be small, 
k < 1 /y/2 leading to type-I superconductivity while keep- 
ing canonical Landau-Ginzburg parameter large, k > 

1/V5 M- 

II. PERSISTENT NON-DISSIPATING 
TOPOLOGICAL CURRENTS AND VORTICES 

Though the idea of non-dissipating topological current 
in vortices was considered long ago [lj] in the context 
of cosmic strings, we are more interested in the recent 
dev elop ments of similar phenomena in high density QCD 
El El El HI and condensed matter systems [IllMl^l- 
The most important result can be formulated as follows. 

Consider QCD with non-vanishing chemical potentials, 
Hl and [Ir, which correspond to two reservoirs of parti- 
cles with different chirality. Due to the chiral anomaly 
a number of interesting macroscopic phenomena could 
occur: induced non-dissipating vector currents on an ax- 
ial vortices, induced axial currents on a vector vortices, 
magnetization of the axial domain wall, induced angular 
momentum by vortex loops, to name just a few. 

Presently we are interested in the phenomena where 
currents are induced in the background of an external 
magnetic field. To be precise, in the chiral limit and zero 
temperature (m q = 0,T = 0), each fermion species q 
makes an additive contribution to the vacuum expecta- 
tion values for the axial and vector currents, 

where $ = J d 2 x±B z (x±) is the total magnetic flux 
through the cross section S, and the fermion current den- 
sities are defined as j A — q^^q, j v = qj 3 q. A few com- 
ments are in order. Formula ^ has a universal nature, 
as it originates from the fundamental quantum anomaly 
and it is not sensitive to whether the magnetic field is 



localized inside of vortex or uniformly distributed over a 
large area S. Also, corrections due to non-zero fermion 
mass and temperature can be explicitly calculated, but 
will not be discussed in the present paper. 

It is worth noting that similar non-dissipating cur- 
rents have been discussed in condensed matter literature 
[3 HI HI- In particular, the expression for anomalous 
supercurrent has been derived for 3 He-A system based on 
the chiral anomaly (see eq.(5.35) and Fig. 23 in ref.[l9(). 
Some suggestions of how these effects can be experimen- 
tally tested were also presented in ref . fl9| . An impor- 
tant feature of the 3 He-A system is the existence of the 
anomalous chiral symmetry, which is not present in 4 He. 
Therefore this phenomenon exists in 3 He-A but not in 
4 He. 

Anomalous supercurrent may also exist in high T c su- 
perconductors with d-wave paring or in a graphene sys- 
tem when the relevant degrees of freedom satisfy the 
massless Dirac equation, with a velocity of order one hun- 
dredth that of light. In these cases, the chiral symmetry 
is obviously present and there is a good chance that non- 
dissipating, topological currents may exist. 

For now we will restrict our discussion to QCD. Note 
that in case of equal chemical potentials, fiR = /i£ = fi, 
the vector current is not induced due to the exact can- 
cellation between left handed and right handed fermions. 
This would be exactly the case in normal nuclear matter 
when 7r meson condensation does not occur or in color 
superconducting phases, such as CFL at asymptotically 
large chemical potentials, where Goldstone modes do not 
condense (see recent review [22] and references therein 
on color superconductivity). 

However, it is known that Goldstone modes are likely 
to condense in nuclear matter [H |24j and will definitely 
condense in color superconductin g p hases for intermedi- 
ate chemical potential [11, [2(|, |27|. (22J. In the cases 
of neutral tt° condensation in nuclear matter \2m and ij 
condensation in the color superconducting phase [271 ] the 
vector current will definitly be induced, as shown by the 
simple argument below. 

The condensation of charged Goldstone mesons 
(tt^, K) in the system is much more complicated, but we 
can expect that the condensation of the charged pseudo- 
scalar Goldstone mesons in the system will lead to dif- 
ferent densities for L and R species (and correspond- 
ing to different effective chemical potentials for R and 
L modes), in which case the vector current will be also 
induced. 

It is not the main goal of this paper to describe all 
possible phases where axial density (and consequently, 
the vector current) can be induced. Rather, we want to 
give a simple argument demonstrating why the condensa- 
tion of a pseudo-scalar Goldstone mode would necessarily 
lead to a difference in densities for left handed and right 
handed species. 

To simplify arguments, we consider a model with just 
a single flavor. Let us assume that we are in a phase 
where the baryon density is non-zero and a neutral r\ 



4 



Goldstone mode is condensed. This implies that our 
ground state can be understood as a coherent superpo- 
sition of an infinitely large number of the Goldstone 77 
mesons. We expect that the ground state of the sys- 
tem is not disturbed by adding one extra 77 meson into 
the system. On the other hand, we can relate the ma- 
trix element with an extra 77 meson to the matrix ele- 
ment without the 77 using the standard PCAC technique, 
(A\0\Atj) ~ i(A\[0,Q 5 ]\A). In the present case the co- 
efficient of proportionality would not be precisely 1/F 
(where F is the Goldstone coupling constant) because 
our Goldstones are in the (77) condensed phase rather 
than in a trivial vacuum. 

Taking \A) to be the ground state and O the baryon 
density operator, one can immediately see that if baryon 
density does not vanish in the ground state, then the axial 
density ([O, Q 5 ]) will not vanish also[51]. This implies 
that densities for left handed and right handed species are 
different and therefore a vector current will be induced. 
We assume this to be the case in what follows. 

We should note that the condensation of the pseudo 
scalar Goldstone mode is not the only mechanism capa- 
ble to produce the asymmetry between R and L modes; 
any P parity violating processes can do the same job. 
The crucial point here is not the ability to produce the 
asymmetry (which is a common phenomenon in neutron 
stars due to the neutrino emission), but the ability of the 
non-dissipating persistent currents to keep this asymme- 
try through the entire volume of the star and deliver it 
to the surface of the star. The phenomenological signifi- 
cance of this is argued in section V. 

As all quarks have non-zero electromagnetic charges, 
once a vector current is induced, an electromagnetic cur- 
rent will also be induced. Our next step is to derive the 
interaction between two superconducting vortices where 
an induced electromagnetic current is present in their 
cores. Our ultimate goal is to understand (at least qual- 
itatively) the changes which will occur in the system due 
to these induced currents. This is the subject of the next 
two sections. 



III. STRUCTURE OF A CURRENT CARRYING 
VORTEX 



glected so E measures free energy per unit length, 



E = d 2 x< — 
I 2m 



■Mbi^wr 



+?I^W| 4 4( Vx A(x)) 2 + ij-A) , (2) 

where ph is the chemical potential of the Cooper pairs 
and a is related to the scattering length I, a = ilTh 1 . The 
tp field describes only the superconducting component of 
the protons. Remember that the fundamental particle 
here is the proton Cooper pair whose mass is actually 
777 = 2m p and its charge is q = 2|e|. The current source 
that was added will be used to model the current flowing 
along the core of the vortex. As discussed earlier, this 
current can be treated as an external electromagnetic 
current. 

The Landau-Ginzburg free energy without the extra 
term for the current source has a symmetry under the 
following gauge transformation, 



A(x) -» A(x)+V^(x) 
t>(x) -> e^ (x) ?/;(x) . 



(3) 



We assume that the external current j is conserved V • j = 
in which case ([3]) obviously remains a symmetry. 

We will choose a form for our current source, j = 
jS 2 (x)z, which models a current at r traveling in the 
z-direction. 

The free energy can be minimized to yield the the equa- 
tions of motion. Minimizing with respect to the vector 
potential A yields a Maxwell equation for our system, 

]- (V 2 A — V (V • A)) = -— jNoether - "j ■ (4) 

An v ' m c 

The right hand side is written in terms of the Noether 
current 

JNoether = J, (V f W> - 7^) - ^A \^\ 2 . (5) 

Minimizing with respect to the order field ip gives 



(6) 



For determining the structure of the vortex we choose 
to place it at the origin and write the ansatz in cylidrical 
coordinates, 



In formulating the problem we assume that the su- 
perconductor is due to non-relativistic proton Cooper 
pairing though similar results are also valid for phases 
where a relatevistic field theory should be used. It does 
not change the qualatative picture described below. We 
start with the two dimensional Landau-Ginzburg free en- 
ergy with term added to model a current source j. As 
discussed earlier, this current is an induced persistent 
electromagnetic current and it couples to the gauge field 
naturally. Dependance along the third direction is ne- 




p(r) e 



i<f> 



A = + z, 



(7) 



(8) 



where riQ = ^ is the density of the superconductor. 
The function p(r) £ (0, 1) describes the profile func- 
tion of superconducting density, p(r) — being no su- 
perconducting material and p(r) — 1 being completely 
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superconducting. If we assume that the current goes 
to zero far from the origin we can use ([5]) to see that 
linv^oo a(r) = 1 and Una r ^oo f{ r ) = 0. The phase of ip 
is chosen to mimic a vortex with winding number n = 1 
and only depends on the (f> coordinate. 

In order to calculate long range interactions between 
vortices we are interested in solutions the equations of 
motion as r — > oo. To decouple our set of differential 
equations it is convenient to define 



p(r) = l + a(r) , 



a (r) = 1 + ra (r) 



(9) 



(10) 



such that a (r) , a (r) — > as r — > oo. Substituting ((TJ) 
and ((5J) into equation Q and linearizing yields the two 
equations, 



where £ = , / ,, h2 is the coherence length. 

In the next section, when calculating the vortex inter- 
actions, it will be useful to "unwind" the phase of the vor- 
tex. This will give a much cleaner solution and is done us- 
ing a gauge transformation ([3]), where <p(x) — — 
The solutions for the field equations become, 



he 



2j 



1 - c a K ( — -r 



(17) 



(18) 



We can now move on to calculating the interactions be- 
tween vortices. 



d 2 a 
dr 2 



Ida 
r dr 



( 4 + )a = 



(11) 



IV. INTERACTION BETWEEN TWO 
CURRENT CARRYING VORTICES 



where A = 



47rq 2 n 



/( r) = ^( X ), 



is the London penetration depth. The 



first equation is the modified Bessel equation of the first 
order. We want a solution that goes to zero as r — > oo so 
we choose the solution to be a modified Bessel function of 
the second kind, a(r) — (j). The second equation 

is just a statement of the Green's function, 



(V 2 - a 2 )K (ar) = -2ir6 2 (x) , 



which implies that f(r) = — -c^Ko (v)- Going back 
through all the substitutions we find that the vector po- 
tential is 



A = 



he 



l 



T 



*■ (D 



C Z A 

C 



G) 



In comparison the standard case without j we see that a 
third term has appeared in the interaction. This attrac- 
tive component will play the crucial role in what follows. 

A similar procedure follows for the solution to the order 
field. Substituting and l|10[) into (O and linearizing 
yields 



Superconducting vortices without currents interact 
(12) through two forces. There is an attractive force caused 
by the superconducting order parameter wanting to have 
one defect instead of two and a repulsive electromagnetic 
force caused by the charges swirling around the vortex. 
Two vortices placed side by side have currents running 
in opposite directions on their nearest sides and it is well 
known that opposite currents repel. 

Suppose there are two wires, placed parallel to each 
other, carrying current. If the currents run in the same 
direction the wires will be attracted to one another. 

(13) 

Now consider that superconducting vortices, instead of 
wires, are carrying the current. There are three forces 
working against each other: the attractive electromag- 
netic force from the current, the repulsive electromag- 
netic force from the gauge field, and the attractive force 
from the order field. If the current were strong enough, 
(14) the attractive force would be strong enough to completely 
cancel the force from the gauge field and the vortices 
would always attract and the superconductor would al- 
ways act like a type-I. 



ld_ 

r dr 



da 
dr 



(15) 



This is a modified Bessel equation of the zeroth order 



which has a solution, <j(x) 
this back we find that, 



C<jK 



1 - c„K n 



V2 



Substituting 



J<t> 



A. Calculations 

The philosophy behind calculating the interaction be- 
tween vortices is to find the energy of the entire system 
and then subtract off the energy of the individual vorti- 
cies as originally outlined in [29] . The technique we will 
use was introduced in [30] and has been used to calculate 
vortex interactions in models with two order parameters 
03 > [SlI- The same philosophy as in [29[ is used but the 
actual calculation becomes much less cumbersome. We 
will reduce the theory to a non-interacting, linear one 
(16) and then model the vorticies as point sources. The inter- 
action energy is then calculated from this linear theory. 
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To make the calculation easier it is useful to use a 
gauge tranformation to remove the phase in t()(x). This 
is described in the previous section and yields the form 
tp = — c). To linearize the theory we expand in p 

and A and keep only quadratic terms to get 



J [a 2m a 



We now add source terms to model the vortices, 



(19) 



£„.,„■<.. = J d 2 x{ra + J-A} , (20) 



where r and are the sources for the fields a and A. 
Minimizing this we get the equations of motion, 



(21) 



X 1 



A = 4nJ . 



(22) 



We want to solve for the sources 3 and r such that a 
and A have the same asymptotic solutions we obtained 
earlier in (17]) and (18]). Using (13]), the derivative of (13]) 
with respect to the radial component r, and the identity 



— Ko(ar) = —aK\(ar) , 
dr 

we can solve for the sources, 



^ 2 Mbr, c2/ \ 
T = 2lT0 (x) , 

m a 



J = 



hcd8 2 {x) - j 



2q dr 



<$>- 



c 



(23) 



(24) 



(25) 



The interaction energy is found by substituting J7" = 
J\ + Ji, A = Ai + A 2 , r = t\ + r 2 and a = <7\ + <r 2 
into the total energy E = E[ ICC + E sourcc and subtracting 
of the energies of the vortices, leaving only cross terms. 
The subscripts 1 and 2 refer to two separate vortices and 
positions xi and x 2 respectively. Using the equations of 
motion we get left over cross terms that are interpreted 
as the interaction energy; 



-Bint = / d 2 x{r 1 a 2 + J 1 ■ A 2 } 



(26) 



Though it is not apparent, the interaction energy is sym- 
metric in the exchange of the subscripts 1 and 2. The 
apparent asymmetry arises when the equations of mo- 
tion for either subscript 1 or 2 are substituted in. Us- 
ing (5]), O, ([2"4"]l and (|2"5)) the interaction energy can be 



i2^ 2 (|x- Xl |)X Vl x ~ x 2l 
ma \ 4 




—6 (x-Xi)i4T 




1 M 2 X 2 JU2 \ „ (d 



(27) 



where d = |xi — x 2 |. To evaluate the second term in the 
integral we made use of equations (|23p and (T5]) . 

If ji and j 2 are set to zero we obtain the interaction 
between gauge vortices without current. The only new 
piece in the interaction is that which comes directly from 
the current. If j\ and j 2 run in the same direction there is 
an attractive force and if they run in opposite directions 
there is a repulsive force. This is the expected result if 
we considered parallel wires carrying current. 

The interaction energy ([27]) determines whether the 
vortices attract or repel and whether we see type-I or 
type-II behaviour in the superconductor. If we set ji = 
]2 = j then there are two cases to explore; one when 
j < and one when j > In the first case the 

first term of ([27]) is positive and we obtain the canon- 
ical behaviour for a superconductor where the Landau- 
Ginzburg parameter decides whether the system exhibits 
type-I or type-II behaviour. 



B. Discussions 

First, let us consider the regime when current is large, 

he 2 



j > 



2q\ ' 



(28) 



and the first term in ([27]) becomes negative. This means 
that there is no longer a repulsive term present in the 
interaction and the vortices will always attract. This 
is the main result of this paper. If the condition (j2"5j) is 
met, all the components of the interaction become attrac- 
tive. While a naive calculation of the Landau-Ginzburg 
parameter suggests it is a type-II superconductor, it ac- 
tually behaves like type-I, and the conflict is resolved. 

In this case an intermediate state will be formed, sug- 
gesting that alternating domains of superconducting and 
normal matter coexist in this regime. This state would 
be the lowest energy state on the phase diagram with 
given B,j satisfying condition (28]). This implies that 
type-II vortices don't form at any time in the neutron 
star's life, even during the short period of cooling when 
the transition to a superconducting state takes place. 
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One should remark that due to the static nature of the 
problem the result ([28]) persists for relativistic systems, 
which may be relevant for study color super conducting 
phases. The equations of motion used to derive (T25|) re- 
main unchanged if we replace \J -^ip — * </>rci and redefine 
the corresponding coupling constants. 

Now let us consider a more realistic case when the cur- 
rent is small, 

In this case we cannot make any precise statements 
within our framework. However, based on experience 
in similar situations in condensed matter systems one 
should expect very dramatic changes to the vortex lat- 
tice when there are currents directed along the external 
magnetic field d3,[3l,[3l,[3i], [11], [13 ■ 

In this literature, the presence of longitudinal currents 
is shown to cause a vortex to develop a spiral- vortex in- 
stability. This instability can be delayed for small cur- 
rents or even stabilized due to the impurities. The les- 
son from these condensed matter systems is that when a 
current aligns with the magnetic field of the vortex the 
properties of the vortex lattice are completely changed 
or destroyed. 

We expect similar behaviour in regions of the neutron 
star where both the Landau- Ginzburg parameter sug- 
gests type-II behaviour and longitudinal currents are in- 
duced. While many features of the system are still to be 
explored, the main massage for the present study that it 
is very likely (similarly to CM studies mentioned above) 
that even small currents ([1]) can completely destroy the 
vortex lattice by replacing it with a new still unknown 
structure. 

It is not our purpose to discuss the rich physics re- 
lated to vortex instabilities resulting from longitudinal 
currents, but rather stress that the resulting state will 
definitely be not the rigid Abrikosov lattice. It is unclear 
what structure will replace the Abrikosov lattice but it is 
reasonable to believe that superconductivity would per- 
sist in this new regime; the energy scales associated with 
currents are much smaller than the superconducting gap. 

It is possible (but not necessary) that the intermedi- 
ate state typical for type-I superconductivity will develop 
and alternating domains of superconducting and normal 
matter would coexist. The size and shape of the domains 
are known to be very sensitive to many things: geometry, 
initial conditions, the method of preparation of a sample, 
boundary conditions, surface effects. As is known, the in- 
termediate state is not in thermodynamic equilibrium in 
the strict thermodynamical sense, but rather depends on 
the history of the system. It is also possible that other 
states, such as Bragg glass phase (3(| would develop, or 
vortex -lattice melting transition would take place [371 ] . 

The exact state is not essential at the moment. What 
is essential is that the Abrikosov lattice is destroyed by 
longitudinal currents. There are many alternative states 



that may replace the Abrikosov lattice. We shall re- 
fer to the absence of the Abrikosov lattice (which is a 
consequence of type-II superconductivity) as a type-I su- 
perconductor that supports the intermediate state even 
though many other phases may result from current in- 
duced vortex instabilities. Therefore, the conflict be- 
tween the precession of a neutron star and the standard 
estimation of the Landau- Ginzburg parameter likely will 
be resolved even when induced currents are small. 



V. CONCLUSION AND SPECULATIONS 

If currents are induced in vortices then we have found 
a mechanism that reconciles the condradiction between 
the precession of neutron stars and the standard pre- 
sumption that there is type-II superconductivity inside a 
neutron star. A sufficiently strong current running along 
the core of the vortex and satisfying inequality ([25]) al- 
lows the vortices to attract even if the Landau-Ginzburg 
parameter indicates they should repel. 

A neutron star would rather form the domain struc- 
ture seen in type-I superconductors rather than the vor- 
tex lattice structure seen in type-II superconductors, thus 
resolving the puzzle. We also argued that even small cur- 
rents along magnetic field can completely change/destroy 
the structure of the Abrikosov lattice. 

A pertinent question is whether these currents can ac- 
tually be induced in neutron stars. The answer depends 
crucially on the details of the specific phase realized in 
the core of a neutron star. As formulated in section III, 
the electromagnetic currents will be induced if the Gold- 
sone modes condense in the presence of a background 
magnetic field. 

If we assume this is the case then there are many ques- 
tions to be considered. How would the magnetic field be 
distributed? What is the fate of these currents? 

A) If the current is large, it is expected that the mag- 
netic field could exist in macroscopically large regions 
where there are alternating domains of superconducting 
(type-I) matter and normal matter - the so called inter- 
mediate state. It has been estimated long ago [9( that it 
takes a very long time to expel a typical magnetic flux 
from the neutron star core. Therefore, if the magnetic 
field existed before the neutron star became a type-I su- 
perconductor (before it sufficiently cooled down), it is 
likely that the magnetic field will remain there. 

The intermediate state is characterized by alternating 
domains of superconducting and normal matter where 
the superconducting domains exhibit the Meissner effect, 
while the normal domains carry the required magnetic 
flux. The pattern of these domains is strongly related 
to the geometry of the problem. The simplest geometry, 
originally considered by Landau [38] , is a laminar struc- 
ture of alternating superconducting and normal layers. 
While precise calculations are required for understand- 
ing of the magnetic structure in this case [53[, one can 
give some simple estimation of the size of the domains 
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using the calculations Landau presented for a different 
geometry. His formula [38| suggests that the typical size 
of a domain is 



(30) 



where R is a typical external size identified with a neu- 
tron star core (R ~ 10 km), while A is the typical width 
of the domain wall separating normal and superconduct- 
ing states. We estimate A ~ A as the largest micro- 
scopical scale of the problem. Numerically, a ~ 10 _1 cm 
which implies that a typical domain can accommodate 
~ 10 4 neutron vortices separated by a distance ~ 10~ 3 
cm. 

B) What is the fate of these currents? It is quite possi- 
ble that the current will travel inside the superconducting 
region in only one direction. The current is conserved, so 
it must make a U-turn and start travelling in the oppo- 
site direction either along the crust of the neutron star 
or through regions of normal matter in the intermedi- 
ate state. This is similar to the case of 3 He-A system 
discussed in ref. [H (see Fig. 23 from that reference). 

C) If the current on the way back travels through the 
normal matter of a different domain (or the crust) then 
large current loops with typical sizes comparable to the 
neutron star radius would be created. These large cur- 
rent loops would induce a coherent toroidal magnetic field 
that, when combined with the poloidal field present in a 
neutron star, would create a non-zero magnetic helicity. 

Such a toroidal magnetic field is apparently neces sary 
to describe the temperature distribution of the crust [3i| • 
It has been also argued long ago that a toroidal compo- 
nent in the magnetic field of a neutron star is necessary 
for stablility of the poloidal magnetic field [ioT |. 

D) The existence of a current making a U-turn near 
the surface of the star may be a key to the understand- 
ing of the long standing problem of neutron star kicks 



, |42|. As is known, pulsars exhibit rapid proper mo- 
tion characterized by a mean birth velocity of 450 ± 90 
km/s. Their velocities range from 100 to 1600 km/s [4l[ 
with about 15% of all pulsars having speeds over 1000 
km/s |13 • Pulsars are born in supernova explosions so 
common theories naturally to look for an explanation in 
the internal dynamics of the supernova. However, three- 
dimensional numerical simulations [H[ show that even 
the most extreme, asymmetric explosions do not produce 
pulsar velocities greater than 200 km/s. Therefore, a dif- 
ferent explanation is needed. 

The origin of these velocities has been the subject of in- 
tense study. Many of the theories involve an assymetry in 
the star's structure, and indeed, many mechanisims are 
capable "in principle" of producing the required asym- 
metry. In the presence of an external magnetic field, the 
neutrinos produced in the star are automatically asym- 
metric with respect to the direction of B. However, the 
most common problem with the suggested mechanisms 
is the difficulty of delivering the produced asymmetry to 
the surface of the star. Only when the assymetry reaches 



the surface of the star may it result in producing the 
proper motion of the entire star. 

Due to their topological nature, the current (fTJ) may be 
capable of delivering the required asymmetry produced in 
the interior of the star to the surface without dissipation. 
Even in a strongly interacting theory, the current (fl]) 
is persistent and non-dissipating. In an environment as 
unfriendly as the dense quark/nuclear matter in neutron 
stars there is still no dissipation due to re-scattering, and 
can be effectively used to deliver information across the 
bulk of the star. 

When the current makes a U-turn on the surface, a 
large amount of momentum (due to photon emission) can 
be transfered to the star. Therefore, this is a unique op- 
portunity to use our topological currents (jTJ) for deliver- 
ing the asymmetry produced in the bulk of the star (e.g. 
due to the Goldstone condensation) to solve the problem 
of neutron star kicks [U, H2| • 

One should notice that the currents may not satisfy 
constraint eq. (|2"8"j) for the explanation of neutron star 
kicks. Indeed, relatively small current is still capable of 
transfering momentum because the U-turn mechanism 
remains operative. 

This asymmetry mechanism is different from most in 
that it does not occur during the supernova, but over 
a long period of time. The momentum transfered from 
each emmited photon is small but if given long enough is 
sufficient to accelerate the neutron star to the observed 
proper velocities. 

E) A different, but likely related phenomena, is the 
recent observation of pulsar je ts [4411 which are apparently 
related to neutron star kicks [45[, |46[ . It has been argued 
that spin axes and proper motion directions of the Crab 
and Vela pulsars are aligned. Such a correlation would 
follow naturally if we suppose that the kick is caused by a 
non-dissipating current, as suggested in D). The current, 
and thus the proper motion, is aligned with the magnetic 
field, which itself is correlated with the axis of rotation. 

As we mentioned above, the U-turn mechanism is nec- 
essarily accompanied by the photon emission (which de- 
livers the momenta required for the neutron star kick). 
It would be very tempting to identify the observed in- 
ner jets [44| with the photons emitted when the current 
makes the U-turn and starts travelling in the opposite 
direction along the crust. In this sense the mechanism 
for the kick is similar to the electromagnetic rocket effect 
suggested previously [46J. 

F) What else could happen with vortices near the 
crust? Reference presents calculations in which vor- 
tices near the crust are grabbed and bundled together 
by the Kelvin-Helmoltz [KH] waves created by the insta- 
bility that arises when there is shear stress between two 
fluids. This bundling of vortices could twist them so they 
no longer line up in an array, but instead form vortex 
loops (vortons). These vortex loops lie in a plane per- 
pendicular to the angular momentum rather than than 
along it. 

Such surface KH instability may explain pulsar 
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glitches [48jj . Spiral vortex instability [32| observed in con- 
densed matter systems may also have some relation to the 
formation of these vortons and to glitches. The helical 
structure of the vortex could expand to the surface and 
transfer its angular momentum to the crust. 

The vortex loops made from superfluid vortices are not 
typically stable (similar to cosmic strings but the 
presence of a current in the core of a superfluid vortex 
and an external magnetic field could make these loops 
stable. If vortons are stable then they could attract to 
each other and form a column. Since they attract very 
close to each other this structure would cease to look like 
a bunch of individual vortices but like a single cylindrical 
vortex sheet which carries a surface current similar to a 
solenoid. Vortex sheets have been studied in the context 
of superfluidity by Landau and Lifshitz [491 ] . 

These and many other consequences of this picture still 



remain to be explored. 
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